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1. Introduction and notation 

The idea that some hidden extra dimensions exist in onr spacetime, dates back to the 
pioneering works of Kalnza and Klein [1]. More recently, the so called ’Braneworld’ 
models have attracted mnch attention in an attempt to explain various theoretical 
difficulties such as the Hierarchy problem PE] and the origin of cosmological inflation 
mu- One of the most interesting of these models from a cosmological point of view, 
is the Randall-Sundrum single brane model (RSII). This model has a single positive 
tension brane embedded in an inhnite AdS bulk. The bulk is assumed to be empty 
except for a cosmological constant and ordinary matter is conhned to the brane. At 
low energies the gravity is localized on the brane via the curvature of the bulk and the 
standard 4D gravity is recovered on the brane, however at high energies the gravity 
’leaks’ into the bulk which leads to the modihcation of the standard general relativity 
|6]. Various aspects and properties of this model has been studied extensively in the 
literature; for example see 

On the other hand, when studying any cosmological theory, a lot of information 
can be obtained by analyzing the behavior of cosmological perturbations and comparing 
it to the current observations of the Cosmic Microwave background and large scale 
structure data. This task has been done in the context of brane cosmology by various 
different authors [9]. As in the standard 4D gravity of general relativity, studying 
the braneworld cosmological perturbations is done by considering two manifolds: A 
background FRW and a perturbed ’physical’ manifold. Perturbation of a quantity 
then dehned by determining the difference between its values in the actual physical 
spacetime and the background reference manifold. The metric perturbations in GR can 
not be uniquely dehned, but depend on the choice of ’gauge’. A gauge essentially can 
be regarded as an identihcation map which corresponds spacetime points in the two 
manifolds. Making a different choice of gauge, may result in the change of the values 
of the perturbation variables. It may also lead to un-physical i. e. not real perturbation 
modes. 

The way around this problem is to work with purely gauge invariant equations. 
This means to hnd combinations of perturbation variables that remain invariant under 
a change of the identihcation map. This method, hrst introduced by Bardeen in 
[To], has the advantage of involving only real unambiguous physical quantities. For 
a comprehensive review of the gauge invariant theory of cosmological perturbations in 
general relativity and its applications in inhation theory and computing the spectrum 
of the cosmic microwave background radiation see references nn [m na E]. 

General relativity, while highly successful in explaining cosmological observations, 
is not the only viable theory of gravity. More importantly some of its problems like the 
singularity problem and the issue of the unihcation of fundamental forces led people 
to try and modify or generalize it. Teleparallel theory of gravity hrst introduced by 
Einstein [15] in an attempt to unify gravity with electromagnetism. This theory in its 
general form is basically the gauge theory for the translation symmetry group [16] . The 
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theory possesses a torsion field which can be regarded as the translational field strength 
corresponding to the coframe field ca. Generally the Lagrangian of this theory lacks 
Local Lorentz invariance. Restoring this symmetry mean imposing some restrictions on 
the form of the Lagrangian which results in a theory which is dynamically equivalent 
to general relativity and is usually called the teleparallel equivalent of general relativity 
(TEGR) in the literature [181 [IH]. Theleparallel gravity and its extensions [201121] has 
generated renewed interest in recent years, in the hope that it may provide solutions 
to some of the cosmological problem like the origin and the nature of the dark energy 
[221123] . It should be stated here that both teleparallel gravity and general relativity 
can be regarded as special cases of a more fundamental gauge theory of gravity called 
Poincare gauge theory which contains both torsion and curvature as translational and 
rotational field strengths respectively [24]. It has been shown that a teleparallel setup 
naturally arises from the low-energy effective field theory induced on D-branes, described 
by a Yang-Mills theory on a fiat noncommutative space which can be regarded as the 
low energy limit of string theory izg. 

In teleparallel gravity one considers a set of n linearly independent vectors e, = eh 
which form a basis in the tangent space on every point of the manifold. The dual of this 
basis d'’ = e^ are the coframes. The dynamical variable in the TEGR theory are 
s, called the tetrads (or pentads in 5D) and they relate anholonomic tangent space 
indices to the coordinate ones. The spacetime metric is not an independent dynamical 
variable here and is related to the tetrad through the relations 




( 1 ) 

( 2 ) 


The inverse of the tetrad is defined by the relation efe^ = Sj. 

Teleparallel geometry, r 4 is defined by the requirement of vanishing curvature. In 
the special case of TEGR the spin connection of the theory is also assumed to be zero. 
This assumption is usually called the absolute parallelism condition and by imposing it 
the connection of the theory will be the Weitzenbock connection defined as [30] 

(3) 

which unlike Livi-civita connection is not symmetric on its second and third indices. 
The curvature of this connection is identically zero and the torsion tensor is 


n. = <id,el - d.el) 


( 4 ) 


Gontorsion tensor which denotes the difference between Livi-civita and Weitzenbock 
connections is 


^ p — 2^'^ P ^ P '^P ' 

and the superpotential tensor is defined as 


S'/" = + 


( 5 ) 


a 


p 06/ 


( 6 ) 
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In correspondence with Ricci scalar, one can define torsion scalar 

T = ^/^TV (7) 

Recently, braneworld models in the context of teleparallel gravity has been studied 
and its setup and fundamental equations has been derived [26l [27j. It is understood 
that the RS setup in TEGR is dynamically different from General Relativity. This fact 
has motivated the authors to study the dynamics of cosmological perturbations in this 
setup. In this paper, by perturbing the coframe components (or pentads in 5D), we 
present the fully gauge invariant cosmological perturbations equations in braneworld 
teleparallel gravity. Unlike the metric perturbation in general relativity which contains 
scalar, vector and tensor modes, pentads can only contain up to vector perturbations. 
The tensor part of the perturbation spectrum in teleparallel gravity originates directly 
from the contorsion tensor, which is the difference between Levi-Givita connection of 
general relativity and Weitzenbock connection of teleparallel gravity. As stated, to make 
the perturbation equations gauge invariant, one must find gauge invariant combinations 
of perturbation variables and rewrite the equations in terms of these parameters. These 
combinations in braneworld teleparallel gravity, in analogy to the Bardeen potentials 
of general relativity, has been derived for both geometric and matter perturbations in 
the following sections. 5D models based on teleparallel gravity are also studied in Refs. 

[2H1129]. 

Except when specifically stated, the notation we use throughout this paper is as 
follows: Lower case Latin letters a, b ,..., i,j ,... run over 1, 2, 3 and label spatial tangent 
space coordinates. The Greek indices a,/3,..., p, z/,.. run over 1,2,3 and refer to the 
spatial spacetime coordinates. The letter / refers to 5D tangent space indices. The rest 
of the upper case Latin Letters M, N, R, L refer to 5D coordinate space indices. 

2. Geometry Perturbations 

So far no trace of fifth dimension has been observed on the brane, therefore one can 
consider it separate from other dimensions. This will give the 5D spacetime, the layered 
structure [24]. Therefore the 5D tangent coordinate system e/ = {ei,e^) can be written 
as 

Cm = C j^ii + e^j^e^ (8) 

Here the first four vectors define the coordinate system on the layer and the fifth 
dimension can be chosen to be normal to the other four dimensions. The above 
orthogonality condition leads us to the point that the expansion of eg cannot contain 
any of 4D e* vectors. Then we have 

6 /. = e\ei + e^^e^ % = e%e5 ( 9 ) 

So the general form of pentad and inverse pentad respectively are 
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and 




K 

0 


K 

hi 


Where 6565 = 1 and + e^^hl = 0 are also satisfied. The capital Latin letters 
indices are 5 dimensional and small Latin and Greek letter run as 0,1,2,3. 

Considering the 5D Anti de sitter background and a 4D FRW brane with fiat spatial 
part. The nonzero unperturbed pentads will be 

eo = n, el = 1, (10) 

Here the Latin and Greek indices are spatial and the overline indices refer to the tangent 
space coordinates. The nonzero unperturbed Weitzenbock connection are 


pO — 
00 ~ 


n 


n 


pQ _ 

i /30--^/3, 


pO — 
05 ~ 


n 


n 


■pa ® ra 

r= -S) 


( 11 ) 


To perturb the geometry of the case, one needs to perturb the pentads. Unlike general 
relativity which has metrics as dynamical variable, here pentads are our dynamical 
variables. By using the orthogonality condition ( 8 ), the general form of the perturbed 
FRW pentad can be written as 


e°o = n(l + A) , = aG^ 

e% = na , = aF^ , 6*5 = 0 


e*o = aR* , e°5 = 0 

4 = 1 + ^ , eV = a(hV + RV) 

( 12 ) 


And the inverse pentads are 

hg° = -(l-A) , hA = --B^ , h,^ = --Ci , hi = -a , h-^ = 0 

n ^ n n 

h-^ = a , = , A5= = 1-^ , hf = ks; - D'‘,) 

a 

(13) 


Where A , (p and a are scalar perturbations and G^, B\ F^, and are vector pertur¬ 
bations and along with scalar ones depend on coordinates (r, x°‘, y). r is the conformal 
time parameter and y is the coordinate of fifth dimension. Note that here we don’t have 
tensorial perturbations since pentad unlike metric is a vector and cannot be perturbed 
tonsorially. The torsion, Contorsion, superpotential and torsion scalar of the above 
pentad is given in appendix A. The calculations has been done up to the first order of 
perturbations. 


To have the gauge invariant geometry, the Lie derivative of the added perturbation 
fields should stay invariant under following coordinate transformations 

^ ^ C = (T,L",y) (14) 

In the language of teleparallel gravity, the Lie derivative of the pentad is 




( 15 ) 
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This important result was first derived in [21] (see also [32]) and is usually called the 
’Stewart Lemma’ in the literature. 

Using the mentioned Stewart lemma, we can calculate the change in the pentad 
perturbation parameters under the gauge transformation as 

A.4 =T+-T (16) 

n 

77 

=-D,T (17) 

a 

AS* =U + -U (18) 

a 

AF, = -D^Y (19) 

a 

AD\ = D,U (20) 

A$ = Y' (21) 

Aa = -Y (22) 

n 

Where a dot denotes differentiation with respect to proper time and a prime denotes 
differentiation with respect to hfth dimension. To consider only the scalar perturbations, 
we dehne the scalar part of the vector perturbations as 

B^ = &B C^ = D^C F^ = D^F = D (23) 

Then the transformations of scalar perturbations are 

• 77 

A + T+-T (24) 

n 

B ^ B + L + -L (25) 

a 

77 

C ^C + -T (26) 

a 

F ^ F + -U (27) 

a 

D^D + L (28) 

$ ^ + W (29) 

1 • 

ot —)■ ot + —Y (30) 

n 

It turns out then that the theory is not gauge invariant. To make it invariant we repeat 
Bardeen’s method. We define new helds out of above old ones in a way that they stay 
invariant under (14). 

Ti = S-i)--F (31) 

a 

= A--C --C (32) 

n n 

\1>3 = $ - a'F - aF' (33) 

^4 = a - -F + -F (34) 

n n 
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These geometric objects, built from the pentad perturbation variables, clearly 
remain invariant under the gauge transformation (16). Comparing with the obtained 
potential in 5D RS setup in general relativity, one can easily hnd the above potentials 
are different from the ones derived in general relativity [9]. 

3. Matter Perturbation 

We repeat the same calculations for the matter perturbations. In RS model of brane¬ 
world, matter is conhned to the brane and bulk can be chosen to contain only a 5D 
cosmological constant. We consider the matter on the brane to be a perfect fluid 
{-p,p,p,p). 

We can generally perturb the energy-momentum as 

5T^ = -p5p, = Tl=p5p5i+pTii (35) 

Where 5p and 5p are scalar, vector and vrf is purely tensor ( = 0). 

Under a 4D coordinate transformation, the perturbed energy-momentum tensor will 
transform like 

5p -^5p + pT + 2T 

Qa Qa + p{La H- La) — pDaT 

a 

6p ^5p + pT--T + lDaL^ 
a 3 

2 5 

'^ajS ^ T ‘^T)(^aLp'j —'JayDgL 

Considering again only the scalar parts, we have qa = Daq , iTap = A^^tt and La 
then 


5p^5p + pT + 2t (40) 

q -)■ q + p{L +-L) - pT (41) 

a 

6p^6p + pT--T+ -AL (42) 

a 3 

71 77 + 2L (43) 

Where (p denotes - and p denotes -. Similar to geometry potentials one can make 
matter potential in a way that they stay invariant under above transformations. The 
gauge invariant matter potentials then are 

n 71 

'^p = 6p-2A--{p-2-)C (44) 

n n 

= q-pB + -pC (45) 

n 

q, Sp--{p--)C-IaD (46) 

n a 6 

^^ = 77-2D (47) 


(36) 

(37) 

(38) 

(39) 
= DaL 
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4. Gauge invariant 5D field equations 

The gravitational action in TEGR is 

Where |e| is the determinant of the pentad e{^. The field eqnations then can be obtained 
by variation of (48) with respect to the pentad 

+= 47rGef(-A5+i(y)S «)(49) 

Where A 5 is the bnlk cosmological constant. Here / refers to 5D tangent space and 
the rest refer to 5D coordinate space and = {—p,p,p,p,0) is the matter which is 
localized on the brane. Varions components of the left-hand side of the pertnrbed 5D 
teleparallel field eqnation (49) is given in appendix B. 

Brane is a 4D hypersnrface which divides the 5D bnlk in two regions. The two sides 
of this hypersnrface is connected via the jnnction conditions. The jnnction conditions 
for the TEGR set-np of RS model has been derived in |26j as 

= 47rG Sj^ (50) 

where um is the nnit vector normal to the brane and can be chosen to be ( 0 , 0 , 0 , 0 , 1 ). 
writing (50) in components, we have 

(J = 0, = 0) : - — (1 - 2 $ - A) - = 47rG-^(l + 6p- 3A) (51) 

an n ^ n-^ 


/ ^ X ^ , 0,' n' 

(/ = 0.iV = «):(- + — 


I I ^ I ^ 

2 n 2 an 


2n 2a 2na^ 


F" + 


1 


n , 

1 




= AnG 


— + -^)R“ + + — 

an n-^ an^ na^ 


F“- -d'^a 

2na 2a? 

‘ gd (52) 


(J = a, A^ = 0 ) : 


1 „ 5a ^ 1 • . n' 5a\^ 1 


n'a 


P 


2^3 2n2^“ 


B„ - 


p pa. ^ 1 

H — 3)^0. - 

n'* an^ 


a^n 


( 53 ) 
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{I = a, N = a) : 


igl-F + _ 6‘d'%) + 

a a > ^ a a n 


i) 

- + (4 + — 

a an 


2d 

a?n 




+ -W” 

^ „o ' _ ) a 




an 


= AncA (^p(i + dp)5- + - pi^D,^ + 2D“ )) (54) 


(/ = 5,A^ = 0) : -idaB^-d^D^ 
n^ 


3® ^ 3a' 

-(1 -2A-^)-\ - a = 0 


an^ 


an 


(55) 


(I = 5,N = a) 




3d 






/3] 


an^ 


B^ + 


4d 

a^n 


—)C'° 
an 




na 




an 


4(9“$ + 4'“ = 0 

a^ a 


(56) 

Note that the left hand sides of the above equations should be evaluated at the position 
of the brane. Using the background field equations and the dehnition of geometric and 
matter potentials (31 — 34) and (44 — 47), we can rewrite various components of the 
field equation (49) in terms of gauge invariant variables. The result for the fully Gauge 
invariant held equations are as follows 


{I = a, N = a) 


a L 


4a“(4$l + ^7^2) - V2(4 Ti + 


2n^a 
h 


9 -^Z/ ''V 9 J-I 9 

an^ an^ 

Aa'n' 

•4d: 


an 


2d 


^iS: - ^$45“ - - ^$25“ 


a^n^ 


a^n 


a^n^ 




a^ 


($; + T')(5:) + -p-At'.d: - ^$',5“ = 0 


a^n 


2 ■" 4’^a 

a^n 


(57) 


(I ^ a, N = 5): 


an' 


d_- 3d_.^ 3 a' ^ ^ ,3a'n' 

2 ^" 4-^ a 2--^ a l + + ^ 


ann 


+ (- 


5a' 


2a3 


2n'a' 

a^n 


3d 


7.2 


2a?n^ 


)daAJA = 0 


)5aTl 

( 58 ) 


(/ = 5 , = a) : - ( 


0- ^ N T /O' ^ N orv T 

— + -^)9 $3 + (7^ + 7^~^)9 ^2 
a^ a^n 2a^ 2a^n 
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. ^ X , fi 3d , • 3d^ ^ 

/ = 5, iV = 5 : -V2(^3 + ^^ 2 ) + (— - = 0 

2 an-^ 2an^ 2a^n^ 


(60) 


/:r Ci^ r. ^ i d'O u'^ ^ o! 

(/ = a, yv = 0) : -dj!, + (— + ^)aa4'4 - ^4-4 


3d „ ^ ad ^ ^ 

+ —da'^i = 0 
2a^n n 


(61) 


(/ = S.AT = 0) : 


Aa'^n 


^ d ^ 1 „9 ^ 4 3a' 2n ', ^, 

(62) 


(/ = 0,1V = 0) : ( 


5d 

4a^n^ 

, an' 


an'’ 


2an^ 

3a'd 

a^n^ 


)V2^2 


2 d ^ , a 

— ^2 + 3( 


an-’ 


a?n 


^ 2a'n ', ^ 

+ - 


an^ 


3n'd 

2an^ 


)^4 


3a'n' 

2na? 




, 1 + 3a' 


na 


)^'i 


= 0 

(63) 


n 


(/ = 0, = a) : + (^ + ^)5“^ 


n 


na- 


a^n 


2an^ 




—a“^' 
2a2 2 


n 


-)a"^2 = 0 


(64) 


^2a2n 2a3' 

The first four of these equations are dynamical equations and the rest act as constraints. 
This set of equations can be solved for four potentials Ti, T 2 , T 3 and 4/4 in the bulk 
subject to boundary conditions provided by the junction conditions (eqs (51-56)). The 
resulting perturbation dynamics is in general different from general relativity. One can 
also bring the junction conditions into gauge invariant form. The results are as follows 


r T AT \ 3a' ^, 6a' ^ 

[I = 0, N = 0) : —T; + —T 3 + 


a'n 


an 


an 


an^ 


T 2 = 47rG—T, 

n^ 


(65) 


2a2 


^ ^“^4+ 


2na 


' 1 

(9aTi = 47rG— 


na-’- 


(I = 0,N = a) : 


( 66 ) 



Braneworld Cosmological Perturbations in Teleparallel Gravity 


11 


(I = a, N = 0) 


1 

2an 


da'^4 + 


3a‘^n' 

2n^ 




1 

an^ 


d ^ 


(67) 


1 2n 1 4 ( 7 ' 

{I = a,N = a): -da^[6: - + -^da ^3 


a 


o?n 


= ATiG—pda{^p + 


( 68 ) 


The other two junction conditions are constraints on geometric potentials and 
background pentad variables at the position of the brane and we use them to bring 
the other four equations into gauge invariant form. These equation act as boundary 
conditions for the bulk perturbation equations. 

5. Conclusion and Discussion 

Gauge properties and symmetries of the teleparallel theory of gravity is essentially dif¬ 
ferent from those of general relativity. As a result, one should be careful when dehning 
perturbations in a manifold with absolute parallelism. The difference in dehnitions of 
perturbations in teleparallel and GR, stems from the fact that when working with the 
Stewart lemma, the Lie derivative is different in two theories. Moreover pentad pertur¬ 
bations in teleparallel gravity can only contain up to vector modes. These items also 
implicate the cosmological perturbations in 4D teleparallel set-up. In 4 dimensions it 
has been shown that the results arisen by the perturbed FRW vierbein are different from 
the ones by perturbed FRW metric [33]. This may bring out interests to study some 
cosmological issues to hnd out more about the conceptual role of torsion in such theories. 

The five dimensional case is more complicated. Because of the different junction 
conditions in this model compared with GR, the results obtained on the brane cannot 
be retrieved by general relativity set-up. In the issue of inflation it has been shown 
that in RS model of TEGR, the inflation index grows faster [2^. Gonsidering these 
points, studying cosmological perturbations of this model could be interesting.In this 
paper, we presented the fully gauge invariant cosmological perturbation equations for 
scalar perturbations in teleparallel gravity by writing down the equations in terms of the 
gauge invariant geometric and matter potentials. These potentials are the teleparallel 
versions of Bardeen’s potentials in hve dimensions. The resulting system of equations 
can be solved to uniquely determine the physical scalar perturbation modes. According 
to the result obtained in this paper, studying the cosmological issues in this setup will 
lead to different conceptual interpretations and also observational predictions. 
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Non-zero components of the torsion 


T% = - [dc^Cg - dyC^ 


7^0 _ I 4/ 

^ 50 - ^ 

= a{d^Fp - dpF^) 

T%a = - da^ 

no = - 
TV = (^ - 

n n n 

no = - aV 

TV = $ - na' 

T'® V Q, /q, 

o' 

rpa 7')^'^ ^ A® 

^ /35 - /3 —^/3 


Contorsion 


(A.l) 


, + -d'nv - 

‘ an^ na^ 

+ 2-nn 


KV = - —+ —5la - -V(l - 2A) - doTV)) 

an ^ n^ V ^ a a >/ 

K'Z = - CbF +(±a + -- 2 - 4 .)« + d'^\ 

a an a a ’ 

‘<"1 = h - -9«o - 4^’« + 4t> + (^ - —+ -C; 

2 L n n"' n^ n n 


K^l = 


‘^ / / 2 

, a n na 




+ 4aoT>'^^' - 


av = - 4 - 9 ^^+ {7^ - —+ 44 )^"'+— 

in^a an 2a^n an an^ an 

K'^i = _ i[(!l + A _ A)b. _ b' 7 + _ h _ A)C-, + !!c-. 

2Lnna a^ a a a 

+ - ^F^ - -nl 
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^ 0 
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T) 

= -{l-2^) + A! 
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K^l 


K^l 


K^l 




1 / 1 


a'’ 


at _ + 4b'- 
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- 1-^F'^ + - 

2 2a^n 

- 4-9^^ + + -F'^ 

$ a' n'a 




, n 
' an^ 


na^ 


:)C' 


2an^ 


pi 




n 




(A.2) 


Superpotential 


Sj^ = 


q O'y _ 


^ = 


‘ ' 2a^ na^ an 

+ - \sid^^Dl. 

" 

_ - (^ + —)SfF^^ - hfF'^^ 

na an^ an a 

- —- Sld,B^) - ^{doD^\) - 5ldoD%) 


an 

2d 
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^(1_2A) + 4-“(?^ + 4-^)« 

an^ n^ n a n n 




a) 
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1 /2a' n' ., 2d Aa' 2n' 


- (-1-h ^4' H- a — (-1- 
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C 50 
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2an 
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Aan^ 2a?n 


5a ^^7 _ ^^7 + ^57 + 
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1 Sr/ 
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2a' 
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Q O 7 
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2an 


2an‘^' 
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'a^n 
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-5 0 
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(A,3) 


Torsion scalar 


T = 


6d2 


6a', a' 


a a 


n . 
n 


4d 


6d , 6d i 

—a H-— 

an an^ 


4a' 2n', _a' a' n\ ^ „a' 


an^ 

+ 


12d^ Qa .2a' 
A-—{— 
a an 

' 


(jPvP 


n 

n^ 


}a 


)D\ + 12-{- + -)<l>-6-A' 
n a a n a 


(A.4) 


Appendix B. 


In this appendix we present the left hand side of the perturbed teleparallel held equations 
in full. 
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- 2aa n^ — a^n $ H-^ )<^a 


d^ 


Tl ^ 

-K9.(A+«.+ |D|) + 5(-^ 


a'n' 


+ TT + 


a^n' 


2a?n^ 


2daB^ - + 35:dgB^ - 2doD^^ + doD^^ - 




f b a? ba? , 5a'a 5a 

rA- —a - 


'2 


\2a3n2 a^n^ 


a^n 


3a' If) an' 2d , 
—D^^ + ^a-^a' + 


2a3 

d 


2a'n' 5a 4a'n'^ 2a' 

+ -^A' 


a^n 


aP 


a^n 


a^ 


a^n^ 


a^n 


a‘^n‘^ 


$)d:- 


20% ^ 2a3 


n 


_n'd ro 

2na^ 


(BA) 
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FqO = dy\--{-% - - - —a5(A + $ + \D\) 

° ^ 2 n zan^ 2a^n 2a?n 2a^n 2an 


® n'/3 


d, -rD 


“ia'a? ?>a’a? 5 a? 

-(1 - 2 $) H- A H-— 

2 ^ ^2 J 2a^n3 


,n' a' 


+ F)d% - (, 


n an ^ ^2a?n^ 2a?n 

^ . , an' 3a'a 3n'a 


3a ^ a ^ 2a'n' 

:^-^)A + 3(^—^- 


- - a„D\) + (^ 


an^ a?n^ 2an‘^' 


r(7W + $) (B.2) 


F! = dy \ - + \{- 

° ^ [ 4an 4 ^ 


1 W _ ^ ^ _ J_r.F , hjii _^ 

4 n^ n^ An^ 4 an^ na"^ 


+ 

a-^n 


^an ^ rr ri-" ^ an*- 

—+ T — 2 ^^^ ~ + *^’+ 1 ^ 1 ) +-^(1 - 2 v 4 ) 

4a^n 4an^ J 2an^ a^n^ 

1 o o ^/9 N 3aa^ ^ ., 3a'a^ 3da' 

- a„Z)%) - —(1 - 2^) + —a + — 


i4)c^ + 

no?’ iaV 


.S, FaF-dmF)- — 

, ^ _2n ' 2n 

,2a' n'. ^ 

“+ n'“+W*“^“ 


AA(1--2A) 

a^n^ 

3a'a^ 3aa^ ^ 
- a H-$ 


a' , 
—a 




(B.3) 


F°. = - €a-,B-<) - ,) 


- -d - 2^ + (- + —)aSt + -a'ifj - -a(4 + A+m 

1 ^ r a? ^ • ,a^n' 5a'a^,^ cA 5da? ^ cA 

+ - 4+ 4 ;;^“ + <i;r + —+ 4 


a!^n', , X „ 1 2 d ^ , 9d^ 

+ 7:72'(““ ) “ “ r,2r,2^[7 a] + 


9d^ 3a'^ lla'^ a'n' 3a'n', 

2an^ a^n ~*~ 4an 4an^ 2n^ 


«' n T^I 

4a2n 4n2 “ 


(B.4) 


Fl = dy\-{% + - -F'^ + , + ^d^A 

I a^ 2an a^ a a^ 2a^ 

^ - i-nF - ^ « - 1 ^ 1 ) 
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+ 


a^n 


do 




5a 


'2 


7a'n' j{' 

2a? 2an 2a 


® /r, n /9 o \ /5a'^ a'n' 

— {dyP — doD ^) + + ::- 1 “ 


A- 
3d 

-2 

an^ 

3d^ 


$) 


3a a 


'^2 


-a 


$ + 


3d^ 

2a? 


an'" 

3 da' 
^an2 


2a^ 2an 2a?n^ 


)$ 


4da', 

a^n 


}a 


3a , 

- a 

an 


(B.5) 


F” = 

a 


-,+^)d^{^ + A+\D\)+dy 

a^ 2na 




2a? 


a 

-t: + 


do 


a^n 


n' 1 , 2a' n' 

5a^d 


a 


a^ 


an 


rv 

2^2 ra 


^ a ■ ,a n 

_F_F —( _ 

A ^ A ^ O, ' 

4n 4n 


a^n 
3a'a? 




4 2 

F®” + t''““ 


4n 

2a' n' 

+ (— H- )d[i3D' 

Zr ry~t ■-( 


a^ 


an 


'2 


an , d ,5a 

3da' ^ dn' ^ 3a'a an' 
4an 2n^ a?n Aan^ 


+ 


fi 

7a'n‘ 


a 


3 4 / 

4 ^ “ 4^“ 4n2 


+ ^(7 _^d A-^ 
2an “ 2a2 “ 4an 


(a - l)5a 




2an Aa?n^ 


3d? , ^ 5d ^ 5d • 

)Fa - -^daa + 


Aa^n 


Aan^ 


a 


(B. 6 ) 


F? = d 


2a3 2a2 


+ 


-d^, 


a-^n 


—d^A - and^^D^ 
2 


a^n' 


dc? a^ • 

+ ( 2 da-^)C'“ + yC'“ + ( 2 a'an+ ^ 


)F" + a^nF'' 


+ 


3da^ 
2 n 
1 


B°‘ — an(9"$ 
2 


a^n 




a 




+ —b''^ 

' An^ An^ An 

“'b“ + (4 


,a?n' 
' 4n 


a'a 


— )C^ _ ^ 

4 4n 4n 4 


a^n^ 


n 

a^n 




4] 


an^ 


:F' 


an 


'4an^ 4a^n^ 


+ 3a'n'^^ 


2na^ 


(B.7) 




^ / t-Ci ^Ct' 


a^n 


a^n 


<90 


d^, 


n‘" 


(a a^ F 


n'a^ 

2 n 


!! 
'4n^ 

3da 


““•(1 _ „)S« + (^ _ 5;^)C» - + aa'-D' 


4n 


4] 


da^ a^ • 

)F“ + ( 2 da - ^)C“ + 

71j Zi Z 


(171 

and^A> - d'^A 

2 
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a^n' 


+ {2a'an H—~ and^'^D^ 


/3] 


+ 


- — 

2na'* an 


, 'n' 

da + 


2a‘^n 


2an^ 


F“ + 


n'a 


2n^a^ 


F° 


n 


5d 5a'^ 2a'n' 

+ 7 ;-1- 


n 

2an 

2 


■C'' 


a 


2n^ 2an^ 2an 


n^ 


d^n 


+ 


n'a' 


+ 


an 


2an^ 2n^ 




Where l-Dl is the determinant of H®,,. 
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